Doping dependence of the vortex-core energy in bilayer films of cuprates 
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The energy needed to create a vortex core is the basic ingredient to address the physics of thermal 
vortex fluctuations in underdoped cuprates. Here we theoretically investigate its role in the occur- 
rence of the Beresinskii-Kosterlitz-Thouless transition in a bilayer film with inhomogeneity. From 
the comparison with recent measurements of the penetration depth in two-unit cell thin films we 
can extract the value of the vortex-core energy p, and show that p scales linearly with T c at low 
doping. 
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One of the most puzzling aspects in the physics of high- 
temperature superconductors (HTSC), that makes them 
substantially different from conventional superconduc- 
tors, is the separation between the fundamental energy 
scales associated to superconductivity: the critical tem- 
perature T c , the zero-temperature superconducting (SC) 
gap A and the superfluid stiffness J s — h 2 p s d±/ / Lm. 1 
Here p s is the superfluid density, measured trough the 
London penetration depth A, and p s d±_ is an effective 
two-dimensional superfluid density, with d±_ a character- 
istic transverse length scale (see discussion below) . While 
in a BCS superconductor the gap formation and the ap- 
pearance of superfluid currents happen simultaneously at 
T c , with A ~ T c , in the HTSC at low doping level the two 
phenomena are essentially decoupled, and T c ~ J s . This 
suggests that the transition can be controlled by phase 
fluctuations, described within an effective XY model for 
the phase degrees of freedom, where J s sets the scale of 
the phase coupling. On a general ground, also the ener- 
getic cost p needed to create the vortex core is connected 
to A~ 2 , i.e to J s . Using standard BCS relations one can 
see that at T = both J s and p are of order of the 
Fermi energy, which is a large energy scale compared to 
T c . Nonetheless, in thin films of conventional supercon- 
ductors J siT) goes to zero as T approaches Tbcs, an d 
p(T) ~ J S (T) is reduced, so that vortex creation becomes 
possible but only at a T c very near to Tbcs- 

In underdoped cuprate superconductors, where the 
BCS picture fails, a clear understanding of the typical 
energy scale which controls the vortex-core formation is 
still lacking. In particular, despite the fact that sev- 
eral experiments 2 ' 3 suggest a predominant role of vortex 
fluctuations, whose occurrence is controlled in a crucial 
way by the value of (U, 5 ~ 8 not much attention has been 
devoted yet to characterize the vortex-core energy p, and 
its relation with T c . In this paper we propose a procedure 
to estimate p using penetration-depth measurements in 
thin films. Indeed, in this case the transition is ulti- 
mately of Beresinkii-Kosterlitz-Thouless (BKT) type, as 
it is signaled by the linear relation between T c and J s , 
and its persistence as doping is changed by electric-field 
effect 9,10 or chemical doping. 11 To clarify the role played 



by the vortex-core energy, let us recall some basic features 
of the BKT physics in a SC film made of few unit cells. 
In Y-based cuprates, whose data we will analyze below, 
there are two strongly-coupled Cu02 layers within each 
cell, which will be considered in what follows as the basic 
2D unit of the layered system, of thickness d = 12 A, cor- 
responding to the unit-cell size in the c direction. Then, 
for a n-cells thick SC film one would expect the system 
to behave as an effective 2D superconductor with areal 
density p 2 s d = ndp s , where p s is the 3D superfluid density 
connected to the penetration depth A. The energy scale 
to be compared to the temperature is then 

j _ h 2 p 2 s d _ ndh 2 p s _ nd <£> 2 
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where $o — (h/2e) is the flux quantum, po the vacuum 
permittivity and we used MKS units as in Ref. [11]. Ac- 
cording to BKT theory, a transition is expected when 
J„/T equals the universal value 2/tt, which gives in terms 
of 1/A 2 (T) the relation: 
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This approach assumes that all the layers of the film are 
so strongly coupled that no mismatch of the SC phase 
in neighboring layers is possible, i.e. the film behaves 
as an effective 2D layer of total thickness nd. How- 
ever, in the case of cuprates it is well known that the 
inter-cell Josephson coupling J± is very weak, so that 
except in a narrow region around T c one expects to 
see the BKT behavior of a single unit cell. As far as 
the superfluid-density behavior is concerned, one would 
then expect that the superfluid density drops already 
when Eq. (2) with n = 1 is satisfied, which is certainly 
the case for totally uncoupled cells. The best situa- 
tion to analyze this effect in real systems is provided 
by the finite-frequency sheet conductivity measurements 
of thin Y 1 . x Ca x Ba 2 Cn 3 7 -s (YBCO) films by Hetel et 
al. . Indeed, in these two-unit-cell thick samples the 
superfluid-density downturn is necessarily between Tgirir 
and Tg^rp, and its exact form depends on the relative 
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strength of J± and vortex-core energy. As we shall see, 
this allows us to extract from the data of Rcf. [11] the 
doping dependence of the vortex-core energy. 

As a starting point we need a model for the BKT SC 
transition in the two-layers system, where for simplicity 
we denote in the following with "layer" a single unit cell 
of the sample. By adopting the formal analogy between 
quantum ID and thermal 2D systems, 6 ' 12 we describe 
each layer (labeled with the subscript 1,2, respectively) 
as a quantum ID sine-Gordon model (J = J dx v s / {2-k)): 



#1,2 = - / C0S(2&, 2 ), 



*» = / 



K{d x e ia f + -(d x <f> h2 f 



(3) 



(4) 



Here the 9i represent the SC phases, and <f>i are the con- 
jugate fields, with [6i(x'), d x (pj(x)] = iirSij5(x' — x), K is 
the Luttingcr-liquid parameter, g u is the strength of the 
sine-Gordon potential, a is the short-distance cut-off and 
v s the velocity of ID fcrmions (which is immaterial in 
the 1D-2D mapping where v s t plays the role of the sec- 
ond spatial dimension, 6,12 t being the imaginary time). 
A vortex configuration for the variable requires that 
§ V0 = ±2ir over a closed loop, i.e. the creation of a 2n 
kink, generated by the exponential of its conjugate field, 
the operator e*^. 12 Thus, the parameter K defines the 
superfluid stiffness and g u the vortex fugacity, as: 



is 7tJ 
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where J = J n =i is the single-layer stiffness. Within the 
standard 2D XY model for the phase the vortex-core en- 
ergy n is controlled by J itself: 5,13 



li XY = 7rJln(2x/2e^)~ yJ, 
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where 7 is the Euler's constant. Even though we will 
treat [i as an independent parameter to be fixed by com- 
parison with the experiments, for the sake of clarity we 
will measure it in multiples of \ixy- The effect of the 
Josephson coupling J± is accounted for by the term: 



H ^ = -^2 Jcoa(e 1 -e 2 ) 7 



(7) 



where g± = wJ±/T. As we shall see below, the inter- 
layer coupling is relevant under renormalization group 
(RG) flow and tends to lock the phases in neighboring 
layers. When this effect dominates over the vortex un- 
binding, the superfluid density is not affected by cross- 
ing the Tgjflp. As the T increases further an additional 
coupling generated under RG flow becomes relevant and 
induces the 2D transition at Tgj^ T . It corresponds to the 
formation of a vortex simultaneously in the two layers: 



g J cos(2(^i+&)). 
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It is then clear that a more convenient basis to study the 
system is given by the symmetric/ antisymmetric fields, 
Osm = (6i ± 02) /V%- The full Hamiltonian then becomes 

H = H° a (K a )+H° s (K s ) + ^ J cos(\/2> s )cos(\/2> a ) 

- ^ J cos(V2fl a ) + ^ J cos(2V20 s ), (9) 

where H° is defined in Eq. (4), K s ^ a = K and the initial 
value of g s is zero, even though it is generated at 0(g 2 ) 
under RG flow (see Eq. (15) below). The superfluid den- 
sity J s is connected to the second-order derivative of the 
free energy with respect to an infinitesimal twist S of the 
phase, d x 0i — ► d x 0i — 6. Since the a field is unchanged 
by this transformation while d x s — > d x s — V25, we im- 
mediately see that J s is given by the asymptotic value of 
K s (£) under RG flow: 



h 2 pf _ K s (£^oo)T 
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The perturbative RG equations for the couplings 
K a , K s , g u , g±,g s can be derived by means of the opera- 
tor product expansion. The result is (see also 6 ' 14-16 ): 
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with I = log(a/oo), where ao, a are the original and RG 
rescaled lattice spacing, respectively. Eqs. (11)-(15) share 
many similarities with the multi-layer case discussed in 
Ref. [6]. If the layers are uncoupled (g± = 0) then 
K s = K a at all scales and the transition occurs when 
g u flows to strong coupling. This happens at K s = 2, 
which corresponds to T^^j, according to the definitions 
(1),(2) and (5). However, when the layers are coupled 
g± grows under RG flow, even if it has initially a small 
value. If the bare couplings -which are T dependent- are 
such that g± becomes of order one before than g u , the 
Josephson coupling term will lock the relative phase a 
in neighboring layers, and g u will flow to zero even if 
K s < 2. However, as soon as K s = 2 the g s coupling 
becomes relevant, signaling the simultaneous vortex for- 
mation in the two layers. This effects makes g u relevant 
as well and the superfluid stiffness K s jumps suddenly 
from the value K s = 1 at Tgj? T to zero. 

As it was discussed in Ref. [6], the range of tempera- 
ture above T^} T where the interlayer coupling allows the 
system to sustain a finite superfluid density is not univer- 
sal, and depends crucially on the value of the vortex-core 
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FIG. 1: (Color online) Temperature dependence of the super- 
fluid density in the inhomogeneous bilayer system, with pa- 
rameter values explained in the text. The downturn is located 
at the intersection with the line n — 1 or n — 2 for uncoupled 
or totally coupled layers, respectively. Inset: expanded view 
for the most underdoped case. The dashed-dotted line is the 
result obtained without averaging over the Jo inhomogeneity. 



FIG. 2: (Color online) (a) 1/A 2 and Howai evaluated at us = 
50 KHz for a single J(T) curve (here /i = 3[ixy)- The finite 
frequency leads to a sharp but continuous decrease of 1/A 2 
across Tbkt, along with a peak in <7i. (b) 1/A 2 and ^ousai 
evaluated at finite frequency using the averaged Jinh- It is 
also shown for comparison the homogeneous curve of panel 
(a) (dashed-dotted line). 



energy, which sets the initial value of the fugacity g u , see 
Eq. (5). Thus, at small fi the system will display a rapid 
downturn already at Tgj^ T , followed by an abrupt jump 
at Tgj^p, while for large values of /x we expect to see only 
the signature of the 2-layers BKT transition at T^ 2 T . 
This is indeed what we observe in the experimental data 
of Ref. [11]. While at higher dopings (25 K< T c < 40 K) 
the superfluid density shows a clear downturn already at 
Trkt-i as the doping is decreased further this signature 
disappears and only the BKT jump at T^ 2 T is visible 
(see inset of Fig. 1). To have a quantitative estimate 
of fi we calculated J s by numerical integration of the 
Eqs. (11)-(15), stopped at a scale I* where g± = s ~ 0(1) 
(we used s = 3), to account for the perturbative character 
of the RG equations. The bare temperature dependence 
of the superfluid density mimics the low-T behavior of 
the data, J(T) = J(T = 0) - aT 2 , with J(T = 0) and 
a extracted from the experimental data well below the 
transition. We also assume that Jj_/ J is independent on 
doping, and we choose J±/J = 10 -3 which is appropri- 
ate in the range of doping considered 17 . Thus, the only 
remaining free parameter is /i, that can be chosen by fit- 
ting the temperature dependence of the data around the 
transition. The result, reported in Fig. 1 as A _2 (T), is in 
excellent agreement with the data of Ref. 11 . 

According to the previous discussion, at Tg^T * ne su ~ 
perfluid density should display a BKT jump, while the 
data of Ref. 11 show clearly a broad tail around the es- 
timated Tbkt- This effect, along with the temperature 
dependence of the real part of the conductivity, cannot 
be attributed only to the finite frequency of the measure- 
ments. Instead, the simplest explanation is the existence 
of a T c inhomogeneity in the sample, which we accounted 
for in the fit of Fig. 1. To clarify this point let us discuss 
for simplicity the pure 2D case for a two-layers thick film. 
In this case, after a transient regime the RG flow of the 



Eqs. (H)-(15) simplify to: 

^~«M, f-P -«.)«., (16) 

with a fixed point at K s = 1, which corresponds to T^ 2 T 
in Eq. (2). The complex conductivity a = o\ + ia% at a 
finite frequency us is given by: 

(7 { UJ ) — ~T2~2 ■ 7 N' 

where s(us) = e 1 + is 2 is complex dielectric constant, due 
to bound and free vortex excitations. Following the dy- 
namical theory of Ambegaokar et al. 18 and Halperin and 
Nelson 19 we estimate these two contributions using the 
RG flow (16) of the BKT coupling, and evaluating e(us) 
at the finite scale — log(r u /a). Here r u = ^/TAD/lj 
is the maximum length probed by the oscillating field, 
where D ~ H/m = 10 16 A 2 /s is the diffusion constant of 
vortices and us the frequency of the measurements. Ac- 
cording to Eq. (17) and (1) A~ 2 = [Lqujoz- Due to the 
finite frequency, the jump of A -2 expected in the us = 
case is replaced by a sharp but continuous drop in a range 
AT^ above T"^j? T . At the same time o\ acquires a fi- 
nite value, with a peak of approximately the same width 
AT W . However, using the value us — 50 KHz correspond- 
ing to the experiments of Ref. 11 the rounding effect on 
A~ 2 and the peak width in at, reported in Fig. 2a, are 
still much smaller than what measured experimentally. 
A more reasonable explanation for the large transient re- 
gion is the sample inhomogeneity. Such inhomogeneity is 
also suggested by tunneling measurements in other fami- 
lies of cuprates, 20 where approximately Gaussian fluctu- 
ations of the local gap value are observed. Even though 
the issue of the microscopic origin of this effect is be- 
yond the scope or this paper, nonetheless we find that an 
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analogous distribution of the superfluid-stiffhess Jo val- 
ues around a given J can account very well for the data 
of Ref. [11]. Thus, we compare with the experiments 
the quantity Ji n h(T) = J dJoP{Jo)J(T, Jo), where each 
J(T,Jq) curve is obtained from the RG equations (11)- 
(15) using a bare superfluid stiffness J = Jq — aT 2 . Each 
initial value Jo has a probability P(Jq) — exp[— (Jo — 
Jo) 2 /2cr 2 ](-\/27rcr) of being realized, where the bare aver- 
age stiffness J(T) = Jo — aT 2 has Jo and a fixed by the 
experimental data at low T, where J exp (T) is practically 
the same as J(T). However, using a variance a = 0.05 Jo 
we obtain a very good agreement with the experiments 
near the transition, as far as both both the tail of A -2 and 
the position and width of <r 1 (w) are concerned, see Fig. 
2b. Observe also that such a variance can be compatible, 
within an intermediate-coupling scheme for the supercon- 
ductivity, with the few times larger distribution of gap 
values (a ~ 0.15A) reported in tunneling experiments. 20 
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FIG. 3: (Color online) Vortex-core energy as a function of T c 
extracted from the fit in Fig. 1. (a) Ratio between fi and the 
fiXY value (6), proportional to the single-layer energy J. (b) 
Absolute value of in K. The dashed line is [s, = 8T C . 

The same finite-frequency analysis is made more in- 
volved in the bilayer case, because the RG equations 
(11)-(15) should be stopped at scales smaller than £ u , 
to prevent the flow of the g± at strong coupling. Thus, 



in Fig. 1 we included only the effect of the T c inhomo- 
geneity by means of the average over P(Jo) discussed 
above (with a = 0.06 Jo for the most underdoped sam- 
ples). This procedure accounts very well for the long tails 
of A~ 2 above the transition (see inset), without affecting 
significantly the estimate of \i. 

Let us now comment on the doping dependence of 
fi(T = 0) reported in Fig. 3. As we said, the measured 
T c crosses over from approximately Tgj^ T at high doping 
to Tgft 2 T at low doping. This is reflected in the dop- 
ing dependence of [if[ixy'- indeed, as we observed for 
the multi-layer case, 6 as fi increases with respect to the 
single-layer stiffness J, the transition moves away from 
Tbkt- -^ ^ s wor th noting that since T c is controlled by 
the competition between the vortex fugacity and the in- 
tcrlayer coupling, in principle the crossover of T c from 
Tbkt t° ^bkt cou ld be obtained also by keeping /i fixed 
and by varying J±. However, to reproduce the data we 
should assume in this case an unlikely increase of Jj_ by 
two order of magnitudes as the doping is decreased, up 
to Jj_/J ~ 10 -1 . Despite the non-universal behavior of 

Hxy, the absolute value of \i reported in Fig. 3b scales 
linearly with T c . This is our central result, which estab- 
lishes a precise relation between the vortex-core energy 
and T c in severely underdoped cuprate superconductors. 

In summary, we analyzed the occurrence of the BKT 
transition in a bilayer system. By means of RG approach 
we computed the temperature dependence of the super- 
fluid stiffness of the bilayer, and we proved the crucial 
role played by the vortex-core energy fj, in controlling the 
transition. Taking into account also the sample inhomo- 
geneity we provided an excellent fit of the experimental 
data of Ref. [11], which allowed us to extract for the first 
time a linear scaling of fi with T c in underdoped YBCO. 
A theoretical understanding of this result is still lack- 
ing and no doubt it would constitute a stringent test of 
microscopic proposals for the underdoped phase. 
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